DESY 01-032 
April 2001 



hep-ph/0104099 



On the Anomalous Dimension of 
the Transversity Distribution h\(x,Q )Q 

(N 

>V 
o3 ■ 



in 



(N 

> 

ON 
ON 

o 
o 



43 

Oh! 

43 



X 



Johannes Bliimlein 



DESY-Zeuthen 
Platanenallee 6, D- 1573 5 Zeuthen, Germany 



Abstract 

We calculate the leading order anomalous dimension of the transversity structure function 
directly using three different methods, the local light-cone expansion in the forward case, 
the non-forward case, and the short-distance expansion of the forward Compton amplitude. 
Our results agree with the original calculation by Artru and Mekhfi [Z. Phys. 45 (1990) 
669], which has been doubted recently. We also comment on the next-to-leading order 
anomalous dimension. 
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1 Introduction 



The transversity distribution hi(x,Q 2 ) QTJ, ^ is one of the central inclusive quantities which 
emerge in the transverse polarized Drell-Yan process and deep inelastic scattering off targets 
with a spin J > 1. This function plays also a crucial rule in proving sum rules for other inclusive 
functions in deep inelastic scattering off polarized targets 0. The transversity distribution is a 
flavor non-singlet function and receives contributions starting with twist-2. The experimental 
measurement of hi(x, Q 2 ) is currently planned at different facilities [§J and the scaling violations 
of this quantity may ultimately be compared with the predictions of perturbative Quantum 
Chromodynamics and thus serve as a novel test of this theory. 

The leading order (LO) splitting function Ph±(z) for the distribution hi(x,Q 2 ) has first been 
calculated by Artru and Mekhfi || using time-ordered or 'old-fashioned' perturbation theory |J 
for the range of momentum fractions z < 1 determining the end-point contribution as in the 
case of the longitudinally unpolarized or polarized flavor non-singlet distribution P NS (z) |7|], 
where the conservation of fermion-number allows this. For the transversity distribution such an 
integral relation does not exist a priori. Therefore one may doubt this procedure, as has been 
done in Ref. [§] recently, where a different result was obtained for the splitting function Ph^z) 
in a direct calculation using a method by Ioffe and Khodjamirian |J. As a result the authors of 
|| conclude that either the splitting function Ph^z) was wrongly calculated in the past or the 
method of Ref. || is inapplicable. 

It is the aim of the present paper to clarify this problem, which is of importance for the correct 
understanding of the evolution of the transversity distribution hi(x,Q 2 ). We will calculate the 
anomalous dimension of the distribution hi(x,Q 2 ) by three different methods. Moreover, it is 
interesting to know, whether one can determine anomalous dimensions using the method ||, 
which may be of special importance for time-like higher order calculations, i.e. in a situations 
where the light-cone expansion is inapplicable. A related question as raised in || on the correct 
end-point behaviour indeed occurs also in next-to-leading order, since the splitting function 
there also was evaluated under some assumptions |TIJ, due to the techniques being applied, 
see 



11, 12 



The paper is organized as follows. In section 2 we recalculate the forward anomalous dimen- 
sion of hi(x,Q 2 ) in leading order using the local light-cone expansion, by calculating the non- 
forward anomalous dimension, from which the forward limit is derived, and the short-distance 
expansion of a related Compton amplitude [EJ. In section 3 we comment on the next-to- leading 
order anomalous dimension of hi(x,Q 2 ) and section 4 contains the conclusions. 



2 LO Anomalous Dimension of hi(x,Q 2 ) 
2.1 Forward- Scattering Anomalous Dimension 

We first calculate the anomalous dimension using the local light cone expansion in the forward 
case, pi = P2 = p. The contributing diagrams are depicted in Figure 1. The Feynman rules for 
the vector-valued operator insertions for two quark and two quark and a gluon line read 

O^.A = a^A u {p.A) n (1) 

n 

0«"\A = ^V^A,^(p'.A) i - 1 A A (p.A)"- i . (2) 

3=1 
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They are denoted by the symbol in Figure 1. Here, p' is the second quark momentum at the 
quark-quark-gluon insertion Figure lb,c, a^ u = (z/2)[7 Ai 7 v — Ti/y^], A p is a light-like vector with 
A. A = 0, t a denotes the SU(3) C generators, g is the strong coupling constant, and A the vector 
index of the additional gluon. 

The anomalous dimension 7 n (5 l (/j)) is obtained by 



d 



(3) 



with 



bare 



(4) 



In leading order the anomalous dimension is a gauge invariant quantity 0. We will work, for 
defmiteness, in the Feynman-gauge. The contributions to the Z-factor in the forward case will 



be denoted by with j 



, e, cf. Figure 1. 




Pi 



P2 



d 



Figure 1: The diagrams for the LO anomalous dimension for (non) forward scattering. (§) are the operator 
insertions Eq. fll],]2]). The curled and straight lines denote the gluon and quark lines, respectively. 

In the limit D = 4 — e — > 4 the term I a vanishes in the Feynman-gauge since the contraction 
of the Dirac-matrices in the numerator results into a term oc e, which cancels the pole of the 
loop-integral. The Z-factor, keeping the quark momentum off-shell, is 



Z n {jj) = 1 + 2a s C F 



In 



— p 




(5) 



with a s = a s j (4ir), a s = g 2 / (47r), and Cp = 4/3. The first and second summand in the brackets 
are due to h + I c and Id + I e , respectively. By Eq. ([$]) the anomalous dimension 



a s C F {l + Aj2~ 



2a s C f 



2SAn) 



dzz n ~ l P hl (z) 



(6) 



is obtained in agreement with the result of Ref. 0, where 5*1(72) = S" =1 (l/j)- Note that we did 
not assume any sum-rule in this derivation. The splitting function (z) reads 



Ph. (z) 



2,a s C f 



-2 + 



(1 



(7) 



'In the MS scheme this holds to all orders. 
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2.2 Non-forward-Scattering Anomalous Dimension 



The diagrams in figure 1 can also be used to calculate the anomalous dimension of hi(x, Q 2 ) in 
the non-forward case, p\ ^ p 2 . The operator insertions corresponding to Eqs. (f[l@) read: 



0> 



ix u o^ e 



-IX.p—K- 



0» x = igt a i u a^x x e^ p+ ^ 



x.k 



(9) 



Here we follow the notation of Ref. ||13|| , where the corresponding insertions for the unpolarized 
scalar operators were given, x is a light-like vector, A the gluon vector index and k the gluon 
momentum. Ki and n 2 are light-cone marks with k± = (k 2 ± fti)/2, and p± = P2 ±pi- We work 
in the Feynman-gauge again and choose = 0, K- — 1 in the following. 

There exist various equivalent representations for the non-forward splitting functions, which 
are known as a—, w—, and near-forward representation, see [13-16]. We will first refer to 
the a-representation, which is directly related to the Feynman-parameter representation of the 
diagrams in Figure 1. Let us denote the different contributions to the non-forward Z— factor, 
cf. Ref. [Jl^ l, by J^, k = a, ...,e. J a vanishes for the same reason as in the forward case. In the 
a-representation the non-forward splitting function is obtained by 



Z7T 



-8(ai) - S(a 2 ) + 5(ai) ( — ) + S(a 2 ) ( — 

\0t2J + \«i 



The [ ] + distribution is defined by 



+ -8{a l )5{a 2 ) . (10) 



dx[f(x,y)]+(f)(x) 



f(x,y)[<f){x) 



(11) 



where / is a distribution out of the space D'[0, 1] x [0,1] and 0(x) e 2)[0,1] a general basic 
function, |fL7ll . The first addend in Eq. (0) is due to Jb + J c and the second due to Jd + J e - 
Eq. ( [3~Q|) agrees with a result in Ref. ||18|| , see also [I5|. There are different possibilities to 
derive the forward anomalous dimension from Eq. (0). One can either use the near-forward 
representation and perform the limit r = x.p-/x.p + — > 0, cf. Appendix D in Ref. |15 



or 



perform a direct integral in the a— or w— representation |14|, 1^], or covert K hl into the local 



representation with respect of the two Mellin-variables n and n', cf. fT3| , [T5[1 and obtain the 
forward anomalous dimension by setting n = n' . 
We first change to the w— representation 



K hl {a u a 2 )Da = K hl { Wl ,w 2 )Dw, 



(12) 



with 



w 2 



«i — a 2 



(13) 



K hl (w u w 2 ) 



—2a s Cp { 5(1 — w 2 + Wi) 



+8(1 —w 2 — wi) 



(1 - w 2 - Wi 
2 

1 - w 2 + w 1 ) + 



-5( Wl )5(l-w 2 )\ , (14) 
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and 

Da = da,\da,2 [9(a.i)9(a 2 )9(l — a\ — a 2 ) + 9(1 — a\)9(l — a 2 )9(ai + a 2 — 1)] (15) 
Dw = —dwidw 2 

x [9(1 +wi- w 2 )9(l -wx- w 2 )9(w 2 ) + 9(1 + w x + w 2 )9(l -w x + w 2 )9(-w 2 )](lQ) 
The forward scattering splitting function is derived by 

P hl (z) = f 1 2 dw 1 K hl (w 1 ,z), (17) 



i-l+z 

cf. Ref. 0, |T§. One obtains from Eq. (0,0) 



P hl (z) = 2a s C F |-2 + j^-j- + - *) j . 

in agreement with Ref. [H. 

Likewise one may first calculate the local non-forward anomalous dimension 7™ 



7 £f = 2a s C F 



2 n' 1 

8nn' -.—9nn' \ , (19) 



n — nn 



where 5 nn > denotes the Kronecker symbol and 9 nn > = 1 for n' < n and otherwise. In the forward 
limit n — n' Eq. (|J) is obtained in agreement with Ref. 0. 

2.3 Derivation of the Anomalous Dimension from a Forward Comp- 
ton Amplitude 

The anomalous dimension for hi(x, Q 2 ) may be also calculated form the short distance expansion 
of the Compton amplitude of an axial-vector and scalar current ji and j 2 directly using the 
method by Ioffe and Khodjamirian ||. This has been tried in Ref. 0] recently. To extract the 
anomalous dimension of the corresponding operator matrix element in this method we start to 
write down the renormalization group equation (RGE) for the Green's function F(q;pi, . . . ,p n ), 
the Fourier transform of (0\T[ji(x)j 2 (0)(f>(xi) . . . 4>(x„ 



[£> + IjM + Ihid) ~ "70?)] F{q; Pl , . . . ,p n ) = , (20) 
with n — 2. The RG-operator is given by 

and •yj k (g), r y(g),'y m (g) are the anomalous dimensions of the currents jk, the outer legs, and 
the mass, respectively. j3(g) denotes the /3— function and \x is the factorization scale. At short 
distances the Green's function has the representation 

F(q; Pl , . . . ,p n ) =J2C k (q)E k ( Pl , . . . ,p n ) , (22) 

k 

with Ck(q) the Wilson coefficients and E k (pi, . . . ,p n ) the hadronic matrix elements. Evaluat- 
ing the Compton amplitude Figure 2 one obtains logarithmically divergent contributions which 
belong to the RGE of the Wilson coefficient, 

P> + lh(9) + liM - loM) Ck(q) = . (23) 
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The corresponding RGE for the hadronic matrix elements E k follows from Eqs. (20,23). Here 
lo k {g) is the anomalous dimension of the composite operator 0^(0) being obtained from 



ji(x)j 2 (0) = J2C k (x)O k (Q) . 



(24) 



7o fe (g) is therefore given by Eqs. (|3]|||6|). The term oc ln(— p 2 / // 2 ) which results from the diagrams 
of Figure 2 




7^75 



P 



P 



d 



Figure 2: QCD 1-loop diagrams to the forward Compton amplitude. 



is 



ic{g) = lh(g) + ij 2 (g) - io k (g) = -2a s C F ■ 2Si(n) . 
The anomalous dimensions of the two currents are PO] 

a s 



(25) 



7 . (a) = -C F a s (s - 1) (s - 3) + 4(s - 15)T F N f + (18s d - 126s 2 + 163s + 291)C A 

I 18 L 

-9(s - 3)(5s 2 - 20s + 1)C F ] } + 0(a 3 s ) , (26) 
where s denotes the spin, i.e. Jj(s=i) = and 



7i(«=o) — — 3Ci?a s + CV 



3 } 6 2. 



a 2 + O(a^) , 



(27) 



which yields the splitting function, Eq. @. This shows that the method of Ref. 0] can be used 
to calculate the anomalous dimensions. 



3 A Remark on the Anomalous Dimension of hi(x,Q 2 ) 
at NLO 

The twist-2 NLO anomalous dimension and splitting function for the distribution h±(x, Q 2 ) was 
calculated in Refs. []IIJ using the cut vertex method of |TT| or the method outlined in flr2"fl . In both 
cases the 5(1 — 2;)-contribution cannot be fixed easily by a direct calculation, see also pi] , since 



the calculation is performed in a physical gauge. In fact, a complete calculation of these terms is 
not yet available. On the other hand, complete calculations of the NLO anomalous dimensions 
as occurring for the various polarized and unpolarized non-singlet and singlet twist-2 parton 
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densities contributing to the structure functions exist in the i?£-gauges, cf. Refs. [22]. Thus these 
direct calculations assure that for the non-singlet '— ' twist-2 quark splitting functions fermion 
number conservation holds, as energy-momentum conservation holds in the singlet sector. 

Knowing this result and working in a class of gauges in which the NLO anomalous dimensions 
are gauge invariant one may take advantage in comparing the 5(1 — z)-terms occurring in the 
quark splitting functions P^ s - (z) and (z). This is now done in the axial gauge. One classifies 
the contributing Schwarz distributions dealing with the H — distribution as an individual one and 
not being expressed due to a composed 'function' of 5— and Heaviside functions in the limit 
e — ► + . Then all the 5(1 — z) terms arise from the self-energy terms, which are independent 
of the operator insertion. Due to this they can be determined without an explicit calculation, 
completing the calculation of the splitting function. The result on the forward scattering splitting 



functions has been extended to the case of non-forward scattering in [23 



4 Summary 

We have shown by three different complete calculations that the calculation of the leading or- 
der anomalous dimension of the transversity distribution hi(x,Q 2 ) as originally being derived 



in Ref. || is correct. Our results also agree with those which have been obtained in Ref. pi 
very recently. It was shown that the method of Ioffe and Khodjamirian can be used to derive 
anomalous dimensions in the present example and similarly for related cases if applied to co- 
efficient functions Ck(x), contrary to the result of |§, which has to be regarded to be wrong. 
The methor of Ref. |24| may thus be of importance for higher order calculations in situations in 
which the light cone expansion does not apply. Furthermore we argued that the 5(1 — z) terms 
of the next-to- leading order calculations, although not yet being calculated directly, are right. 

Acknowledgement. For a conversation I would like to thank W.L. van Neerven and 
D. Robaschik. My thanks are due to D. Boer, B. Pire, and R. Kirschner and W. Vogelsang 
for their interest. 

Note added in proof. After this paper was finished D. Boer pointed out to me that the anomalous 
dimension of h%(x, Q 2 ) was calculated prior to Artru and Mekhfi |5j by Shifman and Vysotsky [p5 | 



and at the same time by Baldraccini et al. p6 . Later also Bukhvostov et al. calculated a Matrix 



element in the quasipartonic approximation which can be interpreted as the LO anomalous 



dimension 27 . I would like to thank R. Kirschner for this remark. We also note that the 



anomalous dimension of hi(x, Q 2 ) at small x was calculated in Ref. 
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